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The quantum mechanical behavior of a particle in a double well defies our intuition based on classical 
reasoning. Not surprisingly, an asymmetry in the double well will restore results more consistent 
with the classical picture. What is surprising, however, is how a very small asymmetry can lead to 
essentially classical behavior. In this paper we use the simplest version of a double well potential 
to demonstrate these statements. We also show how this system accurately maps onto a two-state 
system, which we refer to as a ‘toy model’. 


I. WHERE IS THE PARTICLE? 

Given two wave functions, the one shown in Fig. 1(a) 
and the one shown in Fig. 1(b), which is the correct eigen¬ 
state for a double well potential? 

The actual correct answer is that we do not have 
enough information. We do not have enough informa¬ 
tion because we are only given a picture of the double 
well potential, and not an explicit definition. An ex¬ 
plicit definition of the potential used in Fig. 1 will be 
provided below, and this definition will reveal an asym¬ 
metry not apparent in the figure, that the potential well 
on the right hand side is slightly lower than that on the 
left hand side. The potential is actually drawn with this 
asymmetry, but it is so small as to be invisible to the 
eye on this scale (by many orders of magnitude); the net 
result, however, is that the actual ground state is that 
pictured in Fig. 1(b), with the wave function located en¬ 
tirely in the right side well. Hence, a tiny perturbation 
(the ‘flea’) results in a state very different from the fa¬ 
miliar symmetric superposition of ‘left’ and ‘right’ well 
occupancy (shown in Fig. 1(a)). 


II. THE ASYMMETRIC DOUBLE WELL 
POTENTIAL 

A. Introduction 

The double well potential is often used in quantum me¬ 
chanics to illustrate situations in which more than one 
state is accessible in a system, with a coupling from one 
to the other through tunnelling. [1, 2] For example, in the 
Feynman Lectures the ammonia molecule is used to illus¬ 
trate a physical system that has a double well potential 
for the Nitrogen atom. [3] He used an effective two-state 
model to illustrate these ideas, and in most undergradu¬ 
ate textbooks a two state system is utilized for a similar 
purpose. Here instead we will first focus on a full so¬ 
lution to a double well potential; the features inherent 
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in a two state system will emerge from our calculations. 
Indeed, we will also present a refined two-state model to 
capture the essence of the asymmetry in a microscopic 
double well potential, as we are using in Fig. 1. 
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FIG. 1. (a) State A, “Left” and “Right”, or (b) State B: “Just 
Left” or “Just Right”. In the first state the wave function has 
equal amplitude to be in either well, a linear superposition of 
states so common in quantum mechanics, while in the second 
state the wave function is entirely in the right side well. Which 
is the ground state for the double well potential shown? 

The notion that a slight asymmetry can result in a 
drastic change in the wave function is not new it was 
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first discussed in Refs. [4-6], but in a manner and con¬ 
text inaccessible to undergraduate students. More re¬ 
cently the topic has been revisited [7, 8] to illustrate an 
emerging phenomenon in the semiclassical limit (h —> 0). 
These authors refer to the ‘Flea’ in reference to the very 
minor perturbation in the potential (as in ours) and 
to the ‘Elephant’ [6] (the very deep double well poten¬ 
tial). Schrodinger’s Cat has crept into the discussion 
[7, 8] because the ‘flea’ disrupts the entangled charac¬ 
ter (Fig. 1(a)) of the usual Schrodinger cat-like double 
well wave function (Fig. 1(a)). 

The purpose of this paper is to utilize a very simple 
model of an asymmetric double well potential, solvable 
either analytically or through an application of matrix 
mechanics, [2, 9] to demonstrate the rather potent effect 
of a rather tiny imperfection in the otherwise symmetric 
double well potential. Contrary to the impression one 
might get from the references on this subject, there is 
nothing ‘semiclassical’ about the asymmetry of the wave 
function illustrated in Fig. 1(b). We will show, using a 
slight modification of Feynman’s ammonia example, [3, 8] 
that the important parameter to which the asymmetry 
should be compared is the tunnelling probability; this 
latter parameter can be arbitrarily small. This corre¬ 
spondence applies for excited states as well, along with 
other asymmetric double well shapes. 

B. Square Double Well with Asymmetry 

A variety of symmetric double well potentials was used 
in Ref. 2 to illustrate the universality of the energy split¬ 
ting described there. Here instead we will use perhaps 
the simplest model to exhibit the impact of asymmetry 
- we have checked with other versions and the same 
physics applies universally. This model uses two square 
wells, with left and right wells having base levels Vl and 
Vr, respectively, separated by a barrier of width b and 
height Vo and enclosed within an infinite square well ex¬ 
tending from x = 0 to x = a. For our purposes we assign 
the two wells equal width, w = (a—b)/ 2. Mathematically 
it is described by 

{ oo if a; < 0 or £ > a 
Ho if (a — b)/2 < x < (a + b)/2 
V L if 0 < x < (a - 6)/2 (1) 

V R if (a + b)/2 < x < a, 

and is shown in Fig. 2. We can readily recover the sym¬ 
metric double well by using Vl = Vr. Units of energy 
are those of the ground state for the infinite square well 
of width a, E ° = 7t 2 Ii 2 /(2moa 2 ), where mo is mass of the 
particle, and units of length are expressed in terms of the 
width of the infinite square well, a. We will typically use 
a barrier width b = a/5 so that the individual wells have 
widths w = 2a/5. The height of the barrier then controls 
the degree to which tunnelling from one well to the other 
occurs, and 6 = (Vl — Vr )/2 controls the asymmetry. 


In the following subsection we will work through de¬ 
tailed solutions to this problem, first analytically, and 
then numerically. It is important to realize that these 
solutions retain the full Hilbert space in the problem. A 
‘toy’ model is introduced in a later subsection, and re¬ 
duces this complex problem to a two-state problem. We 
then proceed to illustrate how the two-state problem re¬ 
produces remarkable features of the complex problem, as 
a function of the asymmetry in the two wells. 
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FIG. 2. A schematic of the generic asymmetric square double 
well potential. The well widths are the same but the left- and 
right-side levels can be independently adjusted. 


C. Preliminary Analysis 

Fig. 1 was produced with vo = Vo/E\ = 500, b/a = 
0.2, vl = Vl/Ei = 0, and vr = Vr/E® = 0 in part (a) 
and vr = Vr/E\ = —0.00001 in part (b). The change 
in potential strength compared with the barrier between 
these two cases is 1 part in 50,000,000. This results in 
ground state energies of ei = E\/E\ = 5.827034 and 
ei = Ei/Ei = 5.827025 for the symmetric and asym¬ 
metric case, respectively. Needless to say, either the dif¬ 
ference in potentials or the difference in ground state en¬ 
ergies represent minute changes compared to the tremen¬ 
dous qualitative change in the wave function evident in 
parts (a) vs (b) of Fig. 1. 

The potential used is simple enough that an ‘analyt¬ 
ical’ solution is also possible. The word ‘analytical’ is 
in quotations here because, in reality, the solution to 
the equation for the energy for each state must be ob¬ 
tained graphically, i.e. numerically. While this poses no 
significant difficulty, it is sufficient work that essentially 
all textbooks stop here, and do not examine the wave 
function. [10] 
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Assuming that E < Vq , the analytical solution is where 


Vtr (®) = A sin kx 
ipu(x) =Be KX +Ce~ KX 
ipiii{%) = D sing(a — x) 


k= (2 m{E-V L )/h 2 Y 
re = (2m(V 0 - E)/h 2 Y 
q= (2m(E- V R )lh 2 Y . 


Hnm — 5nm (n 2 E® + (Vl + Vr)w / a + V{)b/ a) 

+ S nm [2V 0 -V L - V R \ sinc(2 n) 

+ (1 - 5 nm ) D nm (V L - Vo + [V R - Fo](-l) n+m ]) 

( 8 ) 

where 


where the regions I, II, and III are depicted in Fig. 2. 
Applying the matching conditions at x = (a± b)/2 leads 
to an equation for the allowed energies: 


D nm = sinc(n — to) — sinc(n + to), 


(9) 


sin (kw) H— cos (kw) ) (sin^ui) + — cos (qw)\ = 

K re J \ re ) 

e~ 2Kb fsin(fcuj) — — cos (kw)j ^sin(gu>) —^ cos (qw)^ , 


( 3 ) 


where w = (a — b )/2 is the width of the individual wells 
and b is the width of the barrier. The matching condi¬ 
tions also provide expressions for the relative amplitude 
D/A: 


D 

c ~Kb 

re sin(fcw) — k cos (kw) 

A 


resin(gic) + qcos(qw) 


or, alternatively, 


D 

p +Kb 

re sin (kw) + k cos (kw) 

A 


Ksin(qw) — qcos(qw) 


( 4 ) 

( 5 ) 


To proceed further, Eq. (3) is solved numerically for the 
allowed energy values. As might be expected the low en¬ 
ergy solutions come in pairs. Once E is known, then so 
are k, q and re, and the D/A ratio can be determined 
through either of Eqs. (4) or (5). In addition, through 
similar relations and normalization, all other coefficients 
in Eq. (2) can be determined. Notice that at first glance 
Eq. (4) suggests that \D\ « |A|, while Eq. (5) suggests 
the opposite. In reality both equations provide the cor¬ 
rect answer, though with limited numerical precision one 
is generally more accurate than the other. Which is more 
accurate depends on whether Vr < Vl or vice-versa. For 
Vr $ Vl, we expect \D\ % \A\. 

Alternatively, we solve the original Schrodinger Equa¬ 
tion numerically right from the start by expanding in 

the infinite square well basis \<f> n (x) = y'^ sin( 2 ^) for 
n = 1,2,3....], and ‘embed’ the double well part in this 
basis. More specifically, we write 


W = 'y]c n \<f> n ), (6) 

n 


and insert this into the Schrodinger Equation to obtain 
the matrix equation, 


E 


H r = E r 

1 x nm'-m J - J n l -n i 


m 


( 7 ) 


sin(7rmu) 

smc(n) =- 

7 xn 


( 10 ) 


As before, w = (a — b)/2 and b are the widths of the 
wells and barrier, respectively. The general procedure is 
provided in Refs. [9] and [2], and the reader is referred to 
these papers for more details. Using either the analytical 
expressions or the numerical diagonalization, the results 
are identical. The advantage of the latter method is that 
the study is not confined to simple well geometries con¬ 
sisting of boxes, and students can easily explore a variety 
of double well potential shapes. 


D. Results 

In Fig. 3 the resulting wave function is shown for a vari¬ 
ety of asymmetries. The result is remarkable. As Vr/E\ 
decreases from 0 to a value of —0.00001, the probability 
density changes from a symmetric profile (equal prob¬ 
ability in left and right wells) to an entirely asymmet¬ 
ric profile (entire probability localized to the right well). 
The other ‘obvious’ energy scales in the problem are the 
barrier height (Vb/E° = 500) and the ground state en¬ 
ergies (Eqs/Ei w 5.827), so these changes in the poten¬ 
tial are minute in comparison. Even more remarkable is 
that as far as the energies are concerned, these minute 
changes give rise to equally minute changes in the ground 
state energy: E GS /E% ~ 5.827034) (« 5.827025) for 
Vr/Ei = 0.0 ( Vr/E j 1 = —0.00001), respectively, while 
the changes in the wave functions are qualitatively spec¬ 
tacular. 


E. Discussion 

That such enormous qualitative changes can result 
from such minute asymmetries in the double well poten¬ 
tial is of course important for experiments in this area, 
where it would be very difficult to control deviations from 
perfect symmetry in a typical double well potential at the 
10 -5 level. Why is this phenomenon not widely dissem¬ 
inated in textbooks? And what precisely controls the 
energy scale for the ‘flea-like’ perturbation that eventu¬ 
ally results in a completely asymmetric wave function 
situated in only one of the two wells [Fig. (lb)]? The 
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FIG. 3. Progression of the wave function as the right well level 
is lowered from 0 (same as the left well level) to —10” 5 (in 
units of Ei = h 2 Tv 2 /(2ma 2 )). When the double well is sym¬ 
metric the probability density (|"0(or) | 2 ) is symmetric (shown 
in red); the degree of asymmetry in the probability density 
increases monotonically as Vr/Ei decreases from 0 to —10 -5 . 
The actual values of Vr used in this plot, in units of E®, are in¬ 
dicated in the legend. Any of these values is absolutely indis¬ 
tinguishable on the energy scale of the barrier ( Vo/E® = 500) 
or the ground state energies ( Egs/E® ~ 5.827). 


answer to the first question is undoubtedly connected to 
the lack of a straightforward analytical demonstration of 
the strong asymmetry in the wave function. Certainly 
Eqs. (4,5) exist, but it is difficult to coax out of either 
of these equations an explicit demonstration of the re¬ 
sulting asymmetry apparent in Fig. (3) as a function of 
lowering (or raising) the level of the potential well on 
the right. To shed more light on this phenomenon, and 
to provide an answer to the second question, we resort 
to a ‘toy model’ slightly modified from the one used by 
Feynman [3] to explain tunnelling in a symmetric double 
well system, and introduced more recently by Landsman 
and Reuvers [7, 8] in a perturbative way. Such a model 
is also used in standard textbooks to discuss ‘fictitious’ 
spins interacting with a magnetic field[ll] and two level 
systems subject to an electric field. [12] 


III. A TOY MODEL FOR THE ASYMMETRIC 
DOUBLE WELL 


An aid towards understanding the results of our ‘mi¬ 
croscopic’ calculations is provided by an ‘effective’ model. 
The tact is to strip the system of its complexity and fo¬ 
cus on the essential ingredients. In this instance, the key 
features amount to whether the particle is in the right 
well, or left well, or a combination thereof. Following 
Feynman, we begin with two isolated wells, each with 
a particular energy level, and with each coupled to the 


other through some matrix element, t: 

Hi/j l = - tij} R 

HijiR = E R tp R - tip L , (11) 

where, in the absence of coupling, the left (right) well 
would have a ground state energy El (E r ), and i/l (ipn) 
represents a wave function localized in the left-side (right- 
side) well. A straightforward solution of this two state 
system results in an energy splitting, as in the symmetric 
case: 


E± 


El + E r 
2 


^ E l — E r 


2 


+ t 2 . 


( 12 ) 


For typical barriers (t << (El+E r )/ 2) and small asym¬ 
metries ( El — E r « (El+E r )/ 2), very little difference 
occurs in the energies, in agreement with the results from 
our more microscopic calculations above. If we define 
S = (El — E r )/2 [= (Vl — V R )/2 for the square dou¬ 
ble well potential], then the ground state wave function 
becomes 


* = 71 f ■ + 7=2 

, (13) 

In the symmetric case we recover the (symmetric) linear 
superposition of the state with the particle in the left 
well, along with the state with the particle in the right 
well (see the remark in Ref. [10]). With increasing asym¬ 
metry, however, say with V R <Vl, i.e. S > 0, the ampli¬ 
tude for the particle being in the right well rises to unity, 
while that for the particle in the left well decreases to 
zero. Our toy model illustrates that the energy scale for 
this cross-over is the tunnelling matrix element, t. This 
energy scale must be clearly present in the microscopic 
model defined in Eq. (1), but it is not there explicitly. 


A. Comparison of the toy model to the microscopic 

model 

To see how well the toy model defined by the two-state 
system in Eq. (11) reproduces properties of the micro¬ 
scopic calculations, we make an attempt to compare the 
results from the two calculations. This is most readily 
accomplished by the following procedure. First, the solid 
curves displayed in Fig. 4 are readily obtained by plotting 
the two amplitudes in Eq. (13), 

^ = 2 ( 1_ 

M 2 = 2 (* + (14) 

as a function of S/t. Then, for one of the results shown 
in Fig. (3), we compute the area under the curve on the 
left; it will correspond to an amplitude |cl | 2 in Fig. 4. 
By placing this value on the appropriate curve in Fig. 4 
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FIG. 4. Plot of the amplitude in each well vs the asymmetry 
parameter, 5/t. The circles indicate the area integrated from 
Fig. 3 on each side of the barrier, as indicated. The large circle 
shows the point used to establish a value of t ss 6.84- 10~ 7 I?i, 
so that the expression in Eq. (14) matches precisely the value 
determined by the more microscopic calculation of the previ¬ 
ous subsection. With this value of t the curves corresponding 
to the expressions in Eq. (14) are also plotted, and the agree¬ 
ment is excellent over the entire range of 5/t. This indicates 
that the ‘toy model’ phenomenology is very accurate. 
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FIG. 5. As in Fig. 3, progression of the wave function as 
the right well level is lowered from 0 (same as the left well 
level) to —3.0 ■ 10 -8 (in units of = h 2 n 2 /(2ma 2 )), but now 
with a barrier height of Vo = 1000E°. When the double well 
is symmetric the probability density (|i/j(x)| 2 ) is symmetric 
(shown in red); the degree of asymmetry in the probability 
density increases monotonically as Vr/E° decreases from 0 
to —3.0 • 10 -8 . The actual values of Vit used in this plot, in 
units of Ei, are indicated in the legend. Even more so than 
before, any of these values is absolutely indistinguishable on 
the energy scale of the barrier (Vo/Ei = 1000) or the ground 
state energies (Egs/Eu « 5.947). 


we are able to extract a value of 5/t and hence an effec¬ 
tive value of t (since 5 = (V^ — Vr)/2 is known). This is 
marked by a large circle in Fig. 4. We have thus identi¬ 
fied a value of t, strictly only defined for the toy model, 
with a specific barrier height and width in the more mi¬ 
croscopic calculations connected with Eqs. (2- 5) or their 
numerical counterparts. We can then vary the value of <5 
(as was done to generate the curves shown in Fig. 3) and 
plot the values of the total probability density in the left 
and right wells as a function of 5/t. The smaller circles 
in Fig. 4 are the results of these calculations, and they 
almost perfectly lie on the curves generated from the toy 
model, thus showing that the asymmetric double well 
system indeed behaves like a two-state system described 
phenomenologically by Eq. (11). We have done this for 
other barrier heights and widths and similar very accu¬ 
rate agreement between the two approaches is achieved. 
We have also carried out such comparisons for excited 
states, and also for other kinds of double wells (e.g. so- 
called Gaussian wells), with similar agreement. 

As an example, in Fig. 5 and Fig. 6 we show results 
analogous to those of Fig. 3 and Fig. 4, but for a dou¬ 
ble well with a barrier with the same width, but with a 
significantly increased height. The sequence of probabil¬ 
ity densities in Fig. 5 is similar to those shown in Fig. 3 
except that the changes in the potential asymmetry are 
orders of magnitude smaller. Fig. 6 then confirms that 


the significantly enhanced sensitivity is due to the signif¬ 
icantly reduced effective ‘hopping’ amplitude t between 
the two wells, such that the asymmetry in probability 
densities as a function of potential asymmetry in Fig. 6 
is as it is in Fig. 4 as a function of 5/t, with both of these 
parameters greatly reduced. 


B. Origin of the coupling t 

The origin of the coupling parameter t in the two-state 
toy model is clearly the possibility of tunnelling that ex¬ 
ists from one well into the other. It is rather involved to 
‘derive’ this parameter t from parameters of the original 
double well potential specified in Eq. (1). In fact it suf¬ 
fices to provide an estimate based only on the symmetric 
case (Vl = Vr = 0), and we provide a brief exposition 
here, following Merzbacher.[13] 

One starts with a variational wave function, of the form 

i>±{x) = ^=(ip L (x) ±xlt R (x)), (15) 

where the ± refers to the ground state (+) or first excited 
state (—), respectively, and the subscript L ( R ) refers to 
the left (right) well, respectively. First taking the two 
wells in isolation, we obtain ipi,(x), for example, with a 
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FIG. 6. As in Fig. 4, plot of the amplitude in each well vs the 
asymmetry parameter, 5/t, but now for a much higher barrier 
potential, Vo/E® = 1000. The circles indicate the area inte¬ 
grated from Fig. 5 on each side of the barrier, as indicated. 
The large circle shows the point used to establish a value of 
t i=3 1.45 • 10“ 9 .Ei, so that the expression in Eq. (14) matches 
precisely the value determined by the more microscopic cal¬ 
culation of the previous subsection. With this value of t the 
curves corresponding to the expressions in Eq. (14) are also 
plotted, and the agreement is excellent over the entire range 
of 5/t. As long as the barrier is sufficiently high to delineate 
two very distinct states (‘left’ and ‘right,’ or ‘dead’ and ‘alive’, 
the two state model works very well. 

solution similar to that in Eq. (2), 

4>l{x) = A sin (kx) 0 < x < w 

= A sin (kw)e-^ x ~ w) x > w, (16) 

and similarly for the well on the right; with all coordi¬ 
nates displaced a distance b to the right, it forms a mirror 
image of the one on the left. This distance can be consid¬ 
ered to be very large at first. In these equations we have 
included an unimportant normalization constant, N± in 
Eq. (15) and A in Eq. (16). The energy splitting between 
the two states is determined by the ‘overlap’ between the 
two wells. A straightforward calculation gives 

tfspiit = J dx ip L (x)Hip R (x) oc e~ Kb « e~ b va> °, (17) 

so as Vq, the height of the barrier, increases for fixed 
width b, the splitting becomes less and less. Note, how¬ 
ever, that the parameter t, first introduced in Eq. (11) is 
proportional to this same quantity: 

t oce~ bVV °. (18) 

This factor is equal to 8.6 x 10 -7 and 2.5 x 10 -9 for 
Vo/E 9 = 500 and 1000, respectively. The actual values 


of t obtained phenomenologically through the fitting pro¬ 
cedure described above are 6.8 x 1CV 7 and 1.5 x 10 -9 , re¬ 
spectively, which tracks very closely these exponentially 
decaying factors. 


IV. SUMMARY 

We have examined the simplest asymmetric double 
well potential and explored the behavior of the wave func¬ 
tion as a function of asymmetry. In the symmetric case 
the ground state is a linear superposition of the particle in 
the left well and the particle in the right well. As the floor 
level of the potential on the right side (' V R ) decreases, 
the probability for the particle to be on the right side 
‘slowly’ increases. The remarkable result of our calcula¬ 
tions is that the energy scale over which the transition 
from symmetric ground state to completely asymmet¬ 
ric ground state can be made arbitrarily small. As our 
‘toy model’ calculation demonstrated, this energy scale 
is controlled by the tunnelling probability between the 
two wells, which is an energy scale that is not obviously 
present in the microscopic parameters (height and width 
of the barrier). In fact, the better well-defined the two 
well system is, the smaller this energy scale. Fig. 3 (or 
Fig. 5) demonstrates this quite dramatically, where im¬ 
perceptibly small asymmetries in the potential give rise 
to a completely asymmetric wave function. There is very 
little indication of this failure from the ground state en¬ 
ergy; instead, it requires a calculation of the wave func¬ 
tion to demonstrate this. 

These calculations serve to demonstrate a number of 
important principles for the novice. First, the numerical 
calculation is ‘simpler’ than the analytical calculation, 
and less likely to lead to error. By this we mean that 
solving for the wave function through Eqs. (3,4,5) is a 
little subtle and, for example, the wrong choice of using 
either Eq. (4) or Eq. (5) can lead to inaccuracies. In 
contrast the numerical solution is straightforward. Then, 
the solution obtained here can be tied to perturbation 
theory. An accurate calculation of the energy can be 
achieved with perturbation theory, but not so with the 
wave function! This ties into the variational principle, 
and teaches the important lesson that a (very) accurate 
estimate of the energy certainly does not imply an even 
qualitatively correct wave function. 
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